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ABSTRACT 



The discovery of source states in the X-ray emission of black-hole binaries and neutron-star low-mass X-ray binaries constituted a 
major step forward in the understanding of the physics of accretion onto compact objects. While there are numerous studies on the 
correlated timing and spectral variability of these systems, very little work has been done on high-mass X-ray binaries, the third major 
type of X-ray binaries. The main goal of this work is to investigate whether Be accreting X-ray pulsars display source states and 
characterise those states through their spectral and timing properties. We have made a systematic study of the power spectra, energy 
spectra and X-ray hardness-intensity diagrams of nine Be/X-ray pulsars. The evolution of the timing and spectral parameters were 
monitored through changes over two orders of magnitude in luminosity. We find that Be/X-ray pulsars trace two dilferent branches 
in the hardness-intensity diagram: the horizontal branch corresponds to a low-intensity state of the source and it is characterised by 
fast colour and spectral changes and high X-ray variability. The diagonal branch is a high-intensity state that emerges when the X-ray 
luminosity exceeds a critical limit. The photon index anticorrelates with X-ray flux in the horizontal branch but correlates with it in 
the diagonal branch. The correlation between QPO frequency and X-ray flux reported in some pulsars is also observed if the peak 
frequency of the broad-band noise that accounts for the aperiodic variability is used. The two branches may reflect two dilferent 
accretion modes, depending on whether the luminosity of the source is above or below a critical value. This critical luminosity is 
mainly determined by the magnetic field strength, hence it dilfers for dilferent sources. 

Key words. X-rays: binaries - stars: neutron - stars: binaries close -stars: emission line. Be 



1. Introduction 

The definition of source states in low-mass X-ray binaries 
(LMXB) and black-hole binaries (BHB) has been a very use- 
ful way to describe the rich phenomenology exhibited by 
these systems in the X-r ay band (Hasinger & van der Klis 1989; 
Muno et all l2002t iGierliriski & Done 2002; Homan & Belloni 



2005t Ivan der Klisl I2006t iKlein-Wolt & van der Kli j 12008 
Kordinget al.ll2008i iDunnet alJ 120101: lBellonill2010 h. A state 
is defined by the appearance of a spectral (e.g., power-law, 
blackbody) or variability component (e.g., Lorentzian) associ- 
ated with a particular and well-defined position of the source 
in the colour-colour (CD) or hardness-intensity (HID) diagram. 
The physical interpretation of these components (i.e., thermal 
and non-thermal Comptonisation, emission from an accretion 
disc, quasi-periodic oscillations, etc.) allows the formulation of 
theoretical models and their confrontation with the observations. 
Different theoretical models differ in the details of the physical 
location and combination of these components. 

As the amount of data provided by X-ray space missions in- 
creased, so did the complexity of the observed phenomenology. 
Although the real picture can be quite complex, two are the ba- 
sic states found in neutron-star and black-hole binaries: a low- 
intensity, spectrally-hard state and a high-intensity, spectrally- 
soft state. The current understanding is that the two states reflect 
the r elative contribution of the accretio n disc on the X-ray emis- 
sion dDone et al ] l2007l:lLin et aPllOOl f). When the disc emission 
dominates, the soft state is seen. Matter in the disc moves in 
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near-Keplerian orbits. The innermost stable orbit of an accretion 
disc around a non-rotating black hole is ~ 37?schw or ~ 100 km 
for a 10 M0 black hole, which is of the same order of magnitude 
as the magnetosphere of a weakly magnetized [B < 10^ G) neu- 
tron star. Therefore, the accretion flows around a low-magnetic 
neutron star and a stellar-m ass black hole in an X-ray binary are 
expected to be very similar (Ivan derKlislll994l) . 

On the other hand, the accretion flow around a strongly mag- 
netized (B ~ 10'^ G) neutron star may be very different because 
in this case the magnetic field begins to dominate the dynamics 
of the accreting flow at much larger distance (r » ^schw) and 
the details of how the magnetic field affects th e accretion flow 
are not fully understood (IBecker & Wo lff 2007). High magnetic 
field neutron stars are found in binary systems orbiting a massive 
(spectral type late O or early B) companion. 

While there have been numerous studies on the character- 
isation of source state s in terms of timing and spectral pa- 
rameters in LMXBs tvan der Klisll2006i, and references therein) 
and BHBs dBeUonil l2oTor and references therein), very little 
work of this kind has been done on high-mass X-ray binaries 
(HMXB), despite the fact that they represent a significant frac- 
tion of the galactic X-ray binary population. Although the dis- 
covery of rapid aperiodic variability (flickering) in the Be X- 
ray pulsar V03324-53 dates back to EXOSAT times dSteUa et all 
1985), timing analysis on accreting pulsars have concentrated 
on the properties and variability of the coherent X-ray emis- 
sion (pulsations) a nd the detection of quasi-periodic oscillations 
(iJames et alJ2010l and references therein). In contrast, the study 
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Table 1. Optical and X-ray information of the systems and the outbursts analysed in this work. The X-r ay luminosity quot ed is 
for the 3-30 keV energy range. The critical luminosity depends on the energy of the cyclotron line (see iBecker et al.ll2012[ and 
Sect. l5.2l i. A "-" sign indicates that either the parameter has not been detected or it is not known. 



Source 


Spectral 


p 

^ spin 


Porb 


e 




QPO 


Distance 


Outburst 


4 


r peak / ^ 


. peak , T 


References 


name 


type 


(s) 


(days) 




(keV) 


(mHz) 


(kpc) 


duration (d) 


(ergs ') 








V 0332+53 


08-9V 


4.4 


34.2 


0.42 


25 


51,223 


7 


100 


3.4 X 10'** 


2.0 


10 


1,2,3 


EXO 2030+375 


B0.5II1-V 


41.8 


46.0 


0.41 


11 


200 


7.1 


155 


1.6 X 10^* 


0.9 


11 


4,5,6 


4U 0115+63 


B0.2V 


3.6 


24.3 


0.34 


12 


2,62 


8.1 


55 


1.4 X 103« 


0.8 


9 


7,8,9 


KS 1947+300 


BOV 


18.7 


40.4 


0.03 




20 


10 


165 


7.1 X 10" 


0.4 


-7+ 


10,11 


lA 0535+262 


BOIII 


105 


111 


0.47 


46 


27-72 


2.4 


45 


6.8 X 10" 


0.4 


<1 


12,13,14 


Swift J 1626.6-5 156 


BOV 


15.4 


133 


0.08 






-10 


-100 


5.2 X 10" 


0.3 




15,16,17 


XTE J0658-073 


09.7V 


160.7 


101 




33 




3.9 


105 


3.6 X 10" 


0.2 


<1 


18,19,20 


lA 1118-616 


09.51I1-V 


405 


24 


0.0 


60 


80 


5 


30 


2.9 X 10" 


0.2 


0.4 


21,22,23,24 


GRO J1008-57 


B1-B2III-V 


93.5 


247.8 


0.68 






5 


50 


2.0 X 10" 


0.1 




25,26,27 



t: Estimated from the break observed in the F - L, diagram, and not from Ery,- (see Sect. l5.2t . 
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[7] Li etal. (2012) 
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[19] Yan etal.. (2012) 
[23] Staubert et al. (2011) 
[27] Naiketal. (2011) 
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Fig. 1. The outburst profiles as seen by RXTE/PCA. Each point represents the average over an observation interval (typically 1000- 
3000 s). The count rate was obtained using the PCU2 in the energy band 4-30 keV. Different symbols denote different source states: 
open circles represent the diagonal branch (DB), while filled circles correspond to the horizontal branch (HB). The flaring activity 
during the peak of the outburst in XTE J0658-073 was analysed separately from the smoother decay. 



of the aperi odic variability and broad-band noise components i s 
very scarce (iBelloni & Hasinge^ll990al:lRevnivtsev et al.ll2009l) . 

Massive X-ray binaries are classified according to the lu- 
minosity class of the optical component into Be/X-ray binaries 
(BeXB), in which the optical companion is a dwarf or subgiant 
star, and supergiant (luminosity class I-II) X-ray binaries. The 



former tend to be transient systems while the latter are persistent 
sources. Be stars are non-supergiant fast-rotating, B-type and lu- 
minosity class III-V stars which at some point of their liv es have 
shown spectral lines in emission dPorter & Riviniusl2003l) . In the 
infrared, they are brighter than their non-emitting counterparts 
of the same spectral type. The line emission and infrared excess 
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Table 2. Journal of the RXTE observations. 



Source 


Proposal 


Time range 


Num. 


Exp. 


name 


ID 


MJD 


obs 


time (ks) 


4U 0115+63 


90089 


53254.1-53280.3 


19 


58.8 




90014 


53282.6-53304.6 


17 


37.6 


KS 1947+300 


50425 


51874.2-51968.6 


33 


82.0 




60402 


51970.2-52078.0 


53 


62.9 


EXO 2030+375 


91089 


53914.9-53998.7 


83 


83.5 




92067 


53999.6-54069.9 


71 


56.6 


V 0332+53 


90089 


53340.3-53365.9 


30 


82.2 




90427 


53367.2-53376.6 


8 


13.1 




90014 


53378.4-53430.5 


42 


79.2 


lA 1118-616 


94412 


54838.3-54838.4 


1 


1.3 




94032 


54841.2-54865.0 


25 


8.9 


XTE J0658-073 


80067 


52932.7-52975.4 


31 


128.0 




504 JU 


3/963. /-33036.4 


43 


/J. 9 


Swift J 1656.6-5 156 


91094 


53723.9-53736.0 


11 


28.9 




91081 


53736.0-53766.6 


18 


68.8 




91082 


53747.7-53835.3 


119 


111.1 




92412 


53840.8-54269.9 


112 


104.2 


lA 0535+262 


94323 


55169.4-55210.5 


27 


126.7 


GRO J 1008-57 


93032 


54426.0-54447.0 


16 


21.2 




93423 


54449.2-54457.2 


4 


6.3 



originate in extended circumstellar envelopes of ionized gas sur- 
rounding the equator of the B star (see Reig 2011, for a recent 
review). 

BeXBs represent one of the most extreme case of X-ray vari- 
ability in HMXBs, displaying changes in X-ray intensity that 
span up to four orders of magnitude. In the X-ray band, BeXBs 
are most of the time in a quiescence state, below the detection 
level of X-ray detectors. When active, the X-ray variability of 
BeXBs is characterised by two type of outbursts: 

- Type I outbursts are regular and (quasi)periodic events, nor- 
mally peaking at or close to periastron passage of the neu- 
tron star. They are short-lived, i.e., tend to cover a relatively 
small fraction of the orbital period (typically 0.2-0.3 Porb)- 
The peak X-ray luminosity during this type of outbursts is 
typically <0.2 the Eddington luminosity for a neutron star 
(LEdd = 1.7xl03« ergs-i) 

- Type II outbursts represent major increases of the X-ray flux. 
The luminosity at the peak of the outburst is close to the 
Eddington limit. These outbursts do not show any preferred 
orbital phase and last for a large fraction of an orbital period 
or even for several orbital periods. The presence of quasi- 
periodic oscillations and the large and steady spin-up rates 
measured during giant outbursts support the formation of an 
accretion disc during type II outbursts. 



The purpose of this work is to investigate the timing and 
spectral variability in BeXBs that display large X-ray intensity 
changes (type II outbursts). Because virtually all BeXBs are pul- 
sars, this work seeks to shed light on the X-ray phenomenology 
of accretion-powered pulsars. 

An attempt to define source states through CD/HID and 
timing analysis was carried out by Reig et al. (2006, here- 
after Paper I) on the hard transient X-ray pulsar V0332+53. 
Subsequently, this kind of analysis was extended to three other 
sourc es: EXO 2030+375, 4U 01 15+63, and KS 1947+300 (lRei3 
12008 1. hereafter Paper II) in an attempt to generalise the results 
of Paper I. These works clearly identified the existence of two 
spectral branches that were called horizontal (HB) and diagonal 
(DB) branches according to the motion of the source in the HID. 
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Fig. 2. Hardness-intensity diagram of KS 1947+300 showing 
two spectral branches: a low-intensity horizontal branch (HB, 
filled circles) and a high-intensity diagonal branch (DB, open 
circles). Arrows mark the flow of time. 



The HB is associated with low-flux and highly variable states, 
whereas the DB appears at high-flux states. 

This paper builds up on the results of Paper II and consti- 
tutes a systematic study of the evolution of the spectral and tim- 
ing variability of all the BeXB detected by RXTE that went into 
type II outburst in the period 1996-2011. Based on the evolu- 
tion and correlations of the spectral and timing parameters and 
the position of the sources in the CD/HID diagrams, we aim to 
define and characterise source states in X-ray pulsars. The main 
difference with respect to Paper II, besides doubling the number 
of systems investigated, is the inclusion of broad-band spectral 
analysis and the finer sampling of the data. In Paper II, the X-ray 
outbursts were divided into several (between 7-9) intervals and 
averaged power spectra for each interval were extracted. Each in- 
terval covered typically a few days. Here, we obtain power spec- 
tra and energy spectra (PCA+HEXTE) for each pointing, which 
is typically few thousand seconds long. This order of magnitude 
increase in sampling resolution allows us to track much faster 
changes, likely to be associated with accretion, whereas the ac- 
quisition of energy spectra allows the search for correlations be- 
tween the spectral and timing parameters. 

Section |2] summarizes the observations and the characteris- 
tics of the instruments used in the analysis. Section [3] describes 
the methodology and techniques used in the data analysis. A de- 
tailed description of the results is given in Sect.|4] Our interpreta- 
tion of the results and the implications of this work are discussed 
in Sect. |5] Finally, conclusions are drawn in Sect. |6] 

2. Observations 

Our aim is to study the characteristics of the population of X-ray 
pulsars with Be companions as a whole. In order to have consis- 
tent and homogeneous data sets, we used data from one single X- 
ray observatory. Thanks to the All-sky Monitor (ASM) and fast 
reaction to unexpected events, the Rossi X-ray Timing Explorer 
(RXTE) has been providing large amount of data of hard X-ray 
transients. Its high time resolution, moderate spectral resolution 
and relatively broad-band response (2-150keV) makes it the per- 
fect observatory to undergo this type of study. 

We analysed archived data obtained by all three instruments 
on board RXTE (Bradt et al. 1993). The data provided by the All 
Sky Monitor (ASM) con sists of daily flux a verages in the en- 
ergy range 1.3-12.1 keV (iLevine et al.lll996l) . The Proportional 
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sc sc sc 

Fig, 3. Hardness (soft colour)-intensity diagram. The soft colour was defined as the ratio 7-10 keV / 4-7 keV. The count rate was 
obtained for the 4-30 keV energy band. Open circles designate points in the DB, while filled circles correspond to the HB. A 
logarithmic scale is used for the count rate (confront with Fig.|2]l. 



Counter Array (PCA) covers the energy range 2-60 keV, and 
consists of five identical coaligned gas-filled proportional units 
giving a total collecting area of 6500 cm^ and prov ides an en- 
ergy resolution of 18% at 6 keV (iJahoda et al.lll996h . The High 
Energy X-ray Timing Experiment (HEXTE) is constituted by 2 
clusters of 4 Nal/Csl scintillation counters, with a total collect- 
ing area of 2 x 800 cm^, sensitive in the 15-250 keV ba nd with 
a nom inal energy resolution of 15% at 60 keV (.Rothschild et alJ 
[T998h . 

Due to RXTE's low-Earth orbit, the observations consist of 
a number of contiguous data intervals or "pointings" (typically 
0.5-1 hr long) interspersed with observational gaps produced by 
Earth occultations of the source and passages of the satellite 
through the South Atlantic Anomaly. Data taken during satel- 
lite slews, passage through the South Atlantic Anomaly, Earth 
occultation, and high voltage breakdown were filtered outQ fol- 
lowing the recommendations of the RXTE tearrQ. 

We have studied all BeXBs that went into outburst (the X- 
ray intensity increased by >100) during the lifetime of RXTE and 
have enough number of observations to allow a meaningful anal- 
ysis. We did not include the BeXB candidate 4U 1901 H-03 be- 
cause no optical counterpart is known for this source and no re - 
Uable estimate of its distance is available (iGallowav et al.ll2()05h . 



' Data were filtered out when the difference between the source posi- 
tion and the pointing of the satellite was greater than 0.02°, the elevation 
angle was smaller than 8° (timing) and 10° (spectra), and contamination 
from electrons trapped in the Earth's magnetosphere or from solar flare 
activity was above 0.1. Good Time Intervals which excluded the times 
of PCA breakdowns were used. 

^ http://heasarc.gsfc.nasa.gov/docs/xte/abc/screening.html 



Because type II outbursts are rare events, there are not many 
studies reporting the evolution of the spectral and/or timing 
parameters over the entire outburst in an accreting neutron 
star. An exception is 4U 0115-1-63, which displayed four giant 
outbursts during the period of time analysed in this work, in 
1999, 2000, 2004, and 2008. The 2000 outburst lacks good 
data coverage. Here we present results for the 2004 event only. 
However, it is worth noticing that the source traced very similar 
HIDs in all outbursts. The two branches, HB and DB, are 
clearly distinguished in correspondence to the X-ray flux level. 
For an analysis of the 2008 outburst the reader is referred to 
iLiet alJ (I2OI2I) and iMuUer et al.1 (120121). The 199 9 and 2004 
events have also been studied bv iNakaiima et alJ (120061) and 
iTsveankov et al ] (I2007h . although these studies focus on the 
dependence of the cyclotron line energy with flux. Spectral 
and/or timing studies covering most of the same giant outbursts 
analysed in this work ha ve also been performed for the follow- 
ing so urces: V0332H-53 (ITsveankov et aL 20 lOt INakaiima et al.l 
l2010h. KS 1947-H30 (Gallowav et alJ l2004. EXO 2030-H375 
(Wilson etal. "2001, lA 1118-616 (NespoH & Reig '201I 
Devasiaetal. 201%, lA 0535-1 -262 (Acciari et al.' 2011]) , 

20061: iNest^oU et al.1 boii . 



XTE J0658-073 
GRO J 1008-57 



(McBrid e_at a 
tNaik et all 



2()llh . 

(Ilcdem et alllloTlllReig et al.ll201 ll) . 



Swift J1656.6-5156 



Figure [T] shows the outburst profiles for each source, cre- 
ated with the count rates of PCU2 in the 4-30 keV band. Table[T] 
summarises some of the properties of the systems and of the out- 
bursts analysed in this work. Table |2]gives the proposal ID, time 
span, the total number of "pointings" or observation intervals, 
and the on-source exposure time of the observations. 
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Table 3. Model photon distributions used in the spectral anal- 
ysis. The spectra were fitted with the function f{E) - A(E) x 
(P{E) + G(E)) X (Cf(E) + Ch(E)) 

Model name as in XSPEC 



A(E) 



PHABS 

CUTOFFPL P(E) = K £ re-£/£c„, 

POWERLAWxHIGHECUT^ P{E) = £-re(£c,„-£)/£roid 



GAUSS 
CYCLABS' 

GABS'" 



G{E) = K-^e 



IE-Ej_) 



Ci(E) 



G,bAE) = e 



(E-Ei)^+wf j 



Cmm{E): 0.03-10 keV interstellar photoelectric absorption 

cross-section. 
A'h: equivalent hydrogen column. 
K: normalization in photons keV"' cm"^ s"' at 1 keV. 
T: power-law photon index, 
iicut : cut-off energy in keV. 
Efa\i : folding energy in keV. 
El : iron line energy in keV. 
cr: iron line width in keV. 
£■, : energy of the cyclotron line in keV. 
£), : depth of the cyclotron line. 
W, : width of the cyclotron line in keV. 
Eq: central energy of the absorption line in keV. 
To: optical depth of the absorption line. 
ctq: width of the absorption line in keV. 
t; Only used in 4U 0115+63. 

i refers to the fundamental (/) or haiTnonic (h). 
f^t; Used in 4U 01 15+63, XTE J0658-073 and 1 A 0535+262 
to fit the "10-keV feature". 



3. Data reduction and analysis 

In this section we describe the different techniques used in the 
data analysis. Heasarc FTOOLS version 6.6. 3 was einploy ed to 
perform data reduction, while XSPEC vl2.6 (lAmaud fl996^ was 
used for spectral analysis. Power spectra were obtained using 
an FFT algorithm implemente d in a Fort r an cod e created by 
us, based on subroutines from Press et al ] (Il996h . Data reduc- 
tion and model fitting were automated so that each observation 
was treated in exactly the same way. The same type of spectral 
and noise components were used for all sources. Therefore, this 
work constitutes the first attempt to perform a consistent, ho- 
mogeneous and systematic study of the variability of accreting 
pulsars during giant outbursts. 



3.1. Colour analysis 

In analogy with optical photometry, one can quantify the broad- 
band X-ray spectral shape by defining X-ray colours. An X- 
ray colour is a hardness ratio between the photon counts in 
two broad bands. X-ray colours were directly obtained from the 
background-subtracted Standard 2 PCA light curves in the fol- 
lowing energy bands: soft color (SC): 7-10 keV / 4-7 keV and 
hard color (HC): 15-30 keV / 10-15 keV. These ratios are ex- 
pected to be insensitive to interstellar absorption effects because 
the hydrogen column densities to the systems, obtained from 
model fits to the X-ray spectra, are in the range 0.6 - 3 x 10^^ 
cm"^. These values only affect significantly the X-ray spectral 
continuum below ~ 2 keV. 



3.2. Spectral analysis 

The spectral analysis was carried out on each individual ob- 
servation using PCA (PCU2 only) Standard 2 mode data and 
Standard (archive) mode from the HEXTE (Cluster A gener- 
ally, and Cluster B for data after December 2005, when Cluster 
A stopped rocking between source and background). Both in- 
struments configurations provide a time resolution of 16s and 
cover their energy ranges with 129 channels. All spectra were 
background-subtracted and dead-time corrected. For each ob- 
servation, the two spectra were fitted simultaneously, covering 
an overall 3-100 keV energy range. A systematic error of 0.6% 
was added in quadrature to the PCA spectra, slightly larger than 
the recommended 0.5% by the instrument tearrQ to obtain more 
conservative estimates of the spectral parameters. No system- 
atic error was added to the HEXTE spectra because uncertain- 
ties are dominated by statistical fluctuations in this instrument. 
In the case of 1 A 0535+262, the system which showed the most 
recent outburst, only PCA data were employed, covering the 
3-60 keV energy range, because cluster B stopped rocking in 
December 2009. The model developed to obtain background 
counts from cluster B to be used with cluster A data is known 
to produce spurious line-like features around 63 keV, which for 
this source lies near a harmonic of a cyclotron resonant scat- 
tering feature. Also, because the signal-to-noise of the HEXTE 
spectra of GRO J 1008-57 was very low, only PCA data were 
used for this source. 

The lack of adequate theoretical continuum models for ac- 
creting neutron stars implies the use of empirical models to de- 
scribe the observations. To compare the results from all sources 
consistently, we used the same (or very similar, if leading to bet- 
ter fit) spectral components. Table |3] gives the photon distribu- 
tion for the different models. To fit the spectral continuum we 
used a m odel composed by a combination of phot oelectric ab- 
sorption (iBalucinska-Church & McCammonlll992l PHABS in 
XSPEC) and a power law with high-energy exponential cutoff 
(CUTOFFPL). In 4U 01 15+63, the CUTOFFPL model left sig- 
nificant residuals in the energy range 10-20 keV. For thi s source 
the POWERLAW x HIGHECUT (see l White etalJl983L and ref- 
erences therein) provided better fits (see Table |3] for the dif- 
ferences between CUTOFFPL and HIGHECUT models). For 
oth er choices of the spectra l continuum the reader is referred 
to Krevkenbohm et al.l(ll999h : ICoburn et al. (2002): Maitra et all 
(12012 , and references therein). In general, the value of the hydro- 
gen column density A^h cannot be well constrained by the PCA, 
whose energy threshold is ~3 keV. When known from previous 
works with higher sensitive instruments below 3 keV, A^h was 
fixed to the known value. To account for uncertainties in the ab- 
solute flux normalisation between PCA and HEXTE we intro- 
duced a multiplicative factor which was fixed to 1 for the PCA 
and let vary freely for HEXTE. 

On top of the spectral continuum a number of discrete com- 
ponents are clearly present in the spectra of accreting pulsars. A 
Gaussian line profile (GAUSS) at 6.4 keV was used to account 
for Fe K fluorescence, whereas cyclotron resonant scattering fea- 
tures (CRSF), which we will refer to as "cyclotron lines ", were 
accounted for with a Lorent zian profile (CYCLABS) (iTanakal 
119861: iMakishima et al.l[T99Q) . The width of the iron line were 
initially let free, but once we checked that there was no signifi- 
cant trend with flux, it was fixed to a value of 0.5 keV. 

The sources 4U 0115+63, XTE J0658-073 and 
lA 0535+262 show significant residuals around 7-10 keV. 



^ http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/pcarmf- 
11.7 
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This feature is known as the "1 0-keV feature", and its origin 
is uncertain (ICobum et al.l 12002). It may appear in emission, 
such as in 4U 0115+63 (Ferrigno et al. 2009; Muller et al] 
120 12h and EXO 2030+375 (iKlochkov et alJ l2007l) and then 
it is modelled with a broad Gaussian emission Une, and it 
is referred to as the "bu mp" model, or in ab s orption, such 
as in XTE J0658-073 (McB ride at afl 120061: iNesnoli et alJ 
l20T2h . In this case, a Gaussian absorption profile (GABS) is 
normally usecfl In EXO 2030+375, this extra broad emission 
component is not needed if a CRSF is added at ~10 keV 
(IKlochkov et alJ 120071: IWilson et alJ l2008l) . In 4U 0115+63, 
the POWERLAWxHIGHECUT model gives comparable 
fits to the combination of the CUTOFFPL model with a 
broad gaussian emission profile. Although the use of the 
POWERLAWxHIGHECUT model has been widely used 
to describe the exponential decay in the spectral contin- 
uum of many accreting pulsars, iri cluding 4U 0115+63 
dTsv^ankov et al.l 120071: iLi et al.l 1201 2h . several studies (see 
e.g. Kreykenbohm et al] 19991) have shown that the cutoff 
energy in HIGHECUT may produce an unphysical break, 
which may be interpreted as a line feature. Because the Ecm 
and the energy of the fundamental cyclotron line occupy the 
same region of the spectrum in 4U 0115+63, the application 
of the POWERLAWxHIGHECUT model to the spectrum of 
this source can lead to eiToneous conclusions regarding the 
variabiHty of the CRSF (Muller etal. 2012). In view of this, 
the spectral analysis of 4U 0115+63 was done using the two 
models. 
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3.3. Timing analysis 

For each observation, a light curve in the energy range 2-20 
keV (PCA channels 0^9) or 2-15 keV (PCA channels 0- 
35), depending on the data mode available, was extracted with 
a time resolution of 2"^ s. The light curve was then divided 
into 128-s segments and a Fast Fourier Transform was com- 
puted for each segment. The final power spectrum is the aver- 
age of all the power spectra obtained for each segment. The fi- 
nal power spectra were logarithmically rebinned in frequency 
and correcte d for dead time effe cts according to the prescrip- 
tion given in iNowak et alJ (Il999l) . Power spectra were normal- 
ized suc h that the integral gives the squared rms fractiona l vari- 
ability ( Belloni & Hasinger' ll990bl: [Miyamoto et alj[T991ft . To 
have a unified phenomenological description of the timing fea- 
tures within a source and across different sources, we fitted the 
noise components with a function consisting of one or multiple 
Lorentzians, each denoted as L,-, where / determines the number 
of the component. The characteristic frequency v^ax of L, was 
denoted as v, . This is the frequency where the component con- 
tributes most of i ts variance per loga rithmic frequency interval 
and is defined as dBelloni et alj|2002h 



^(v^ + (FWHMI2Y) = vo Vl + 1/42' 



(1) 



where vq is the centroid frequency and FWHM is the Lorentzian 
full-width at half maximum. Note that Vmax 5^ vq. The broader 
the noise component, the larger the difference between Vmax and 
vq. The quality factor Q is defined ?& Q - vq/FWHM, and is 
used as a measure of the coherence of the variability feature. 
Peaked noise is considered to be a quasi-periodic oscillation 
(QPO) ife>2. 



Note that the GABS model was initially proposed to fit CRSFs 
dSoong et alJl99(]l : ICobum et al.ll2002h 



Fig. 4. Relationship between the X-ray colours (SC= 7-10 
keV/4-7keV, HC=15-30keV/10-15keV) and photon index. The 
SC and the HC of the horizontal branch decrease as the spectrum 
becomes softer, as expected. In contrast, a positive correlation is 
seen in the HC of the diagonal branch. 



The peaks from the neutron star pulsations were fitted to 
Lorentzian functions with the frequency fixed at the expected 
value and width equal to 0.001 Hz (approximately, the inverse 
of the timing resolution of the power spectra). The results from 
the use of Lorentzians are best visualized using the vxPy repre- 
sentation, where each power is multiplied by the corresponding 
frequency. 



4. Results 

In this section we present the results of our analysis. First, we 
performed an independent analysis with each one of the three 
techniques described above, namely, colour-intensity diagrams, 
energy spectra and power spectra. We searched for correlations 
of the colour, spectral and timing parameters as a function of 
X-ray flux and between various spectral parameters. Then, we 
investigated whether it is possible to establish relationship be- 
tween parameters across different techniques. The parameters 
involved in these correlations are the soft (SC) and hard (HC) 
colours, the photon index (F), the cutoff energy (ficut), the en- 
ergy (Epf.) and intensity (/pe) of the iron line and the character- 
istic frequency (v,) and fractional amplitude of variabihty (rmsi) 
of the broad-band noise components. 

To quantify the strength of the relationships between differ- 
ent parameters we used the Pearson's correlation coefficient r, 
whereas the probability p was derived to assess the significance 
of the correlation, p expresses the probability that the observed 
correlation happened by chance. The smaller the p, the more 
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Fig. 6. Power-law photon index as a function of X-ray luminosity. Open circles represent the DB (supercritical regime), filled circles 
the HB (subcritical regime). Orange data points in XTE J0658-073 correspond to the flaring episode during the peak of the outburst 
(see Fig. [TJ. In 4U 01 15+63, red squares represent the results from the "bump" model (see text). The Eddington luminosity for a 
neutron star was taken to be 1.7 x lO-'*' erg s"'. 
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Fig. 5. Average energy spectra of the horizontal and diagonal 
branches at various flux levels. Black and red symbols represent 
PCA and HEXTE data, respectively. The quoted X-ray luminos- 
ity corresponds to the 2-100 keV range. 



significant the relationship. If p < a, the probability that the 
relationship happened by chance is small. (1 - or) x 100 is the 



confidence level (typically 95% or 99%). p was computed from 
the t - value defined as 



N-2 
1 -r2 



(2) 



where N is the number of data points. N-2 gives the number of 
degrees of freedom and p - T(N-2, t), where T is the Student's 
t distribution (two-tailed). 

Table |4] presents a compilation of the Pearson's correlation 
coeflicients and their significance for most of the samples dis- 
cussed in this paper. To have an approximate perception of the 
quality of the fits. Columns 7 and 8 in this table give the mean 
and standard deviation of the reduced values obtained from 
the fits to the energy and power spectra, respectively. In a few 
cases, most notably in lA 0535-1-262, the best spectral fits gave 
reduced xl slightly lower than 1 . In those cases, however, the 
residuals confirmed that all the components employed by the 
model were necessary to describe the data. The exclusion of any 
of them, resulted in an overall not acceptable fit due to marked 
residuals that could be easily fitted with known components. For 
example, the removal of the cyclotron line in the higher-flux 
spectra of lA 0535+262 increased the value of;^'^ from 50-70 
for 77 degrees of freedom to above 200 for 8 1 degrees of free- 
dom. Likewise, removing the broad Gaussian component reulted 
m an mcrease of;t-^ to 90-100 for 80 de grees of freedom. 
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Fig. 7. Cutoff energy as a function of the photon index. Open circles represent the DB, filled circles the HB. Orange data points in 
XTE J0658-073 coiTespond to the peak of the outburst. In 4U 01 15+63, red squares represent the results from the "bump" model. 



10'' 




Fig. 8. Flux of the iron line as a function of the X-ray continuum flux. We assume Ledd = 1.7 x 10^^ erg s Open circles represent 
the DB and filled circles the HB. 
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Table 4. Pearson's correlation coefficients for the vaiious samples. All correlations are significant at > 99%, unless stated otherwise. 
The Lx - v\ and Lx - he relationships were fitted with a power law function (linear correlation in log - log scale). A " — " indicates 
that no data points exist. Also shown (columns 7 and 8) is the mean and standard deviation of the reduced of the fits. In 
4U 0115+63, quoted value corresponds to the POWERxHIGHECUT model. Using the "bump" model the;^^ reduces to 1.0 ± 0.2 
for this source. 
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4. 1 . Hardness-intensity diagrams 

Paper I and Paper II showed that hard X-ray transients exhibit 
two distinct spectral branches in the HID that were called the 
horizontal branch (HB) and the diagonal branch (DB) in Paper II. 
As in LMXBs, these names were adopted because of the pattern 
that the source traces in the HID. Figure|2]shows a characteristic 
example of HID, that of the source KS 1947+300, where the two 
branches can be clearly distinguished. The average intensity in 
the HB is always lower than that of the DB. 

The HB appears horizontal in Fig. |2] because a linear scale 
was used for the Y-axis. Had we used a logarithmic scale, the 
HB would also have a diagonal pattern, although with opposite 
slope (Fig. [3]). In this work we shall use a logaiithmic scale be- 
cause it facilitates the comparison i) with other fainter accreting 
X-ray pulsars, as the behaviour of the points that populate the 
HB stands out clearer and ii) with the HID of BHB and LMXB, 
as for this type of binaiies a logarithmic scale is the usual way 
to represent the count rate. Nevertheless, for the sake of con- 
sistency with previous work, we shall keep the term horizontal 
branch to designate the low-intensity state and diagonal branch 
for the high-intensity state. 

The hardness-intensity diagrams (HID) of all the sources 
analysed in this work, where the count rate in the 4-30 ke V band 
is plotted as a function of the soft colour, are shown in Fig. |3] 
An important result that emerges from the colour analysis is that 
only the brightest sources display the DB (columns 10 and 11 
in Table [T] give an estimate of the peak luminosity). However, 
the X-ray luminosity does not seem to be the only parame- 
ter tiiggering the transition between branches. KS 1947+300 
and lA 0535+262 exhibit similar peak luminosity but only 
KS 1947+300 displays the two branches. We shall come back 
to this different behaviour in Sect. |5] and show that the magnetic 
field strength plays a crucial role in determining the critical lu- 
minosity at which the transition takes place. 

Due to their unpredictable nature, the initial stages of the 
rise of the outbursts lack a proper coverage. Normally, non- 
scheduled observations are triggered once the source is relatively 
bright and the reaction timescales are often similar to the time 
scale of the rise. Only one system, KS 1947+300, was observed 
for the complete duration of the outburst. In all other systems, 
the first data point is located, in the best cases, half way through 
the rise. We shall take KS 1947+300 as an example of the mo- 
tion of the source in the HID. As the X-ray flux increases the 



SC increases, i.e., the source moves right and depicts the HB 
(Fig. |2] see also Fig. O. When the luminosity reaches a critical 
value, the source enters the DB by making a sudden turn in the 
HID: it moves left and the SC starts to decrease. The peak of 
the outburst corresponds to the softest state of the DB. Then, as 
the flux decreases, the source moves back following the same 
track but in the opposite direction. Below the critical luminosity, 
the source enters the HB and the SC begins to decrease again. 
Arrows in Fig. |2] indicate the flow of time. In the sources that 
exhibit only the HB, the peak of the outburst corresponds to the 
hardest spectiTim (largest SC). 

Unlike the X-ray spectra of BHBs and LMXBs where the 
hard emission (i.e., above 10 keV) is well characterised by a 
simple power law component, the higher energy part of the X- 
ray spectra of BeXB is affected by extra components such as ex- 
ponential cutoffs and cyclotron lines and their harmonics. As a 
result, the X-ray colours do not always represent a reliable mea- 
surement of the spectral slope. The HC is expected to be more 
affected by the distortion of the continuum because these extra 
components mainly appear above 10 keV. If the spectrum could 
be well represented by a single power law, then we would expect 
the X-ray colours to decrease as the emission becomes softer, 
that is, we would expect an anticorrelation between the photon 
index and the colours. The SC in both branches and the HC in the 
horizontal branch display this anticorrelation, albeit with large 
scattering in the case of the HC (Fig.|4]i. However, the evolution 
of the HC in the DB is the opposite to what it is expected, that is, 
the HC increases but the gamma also increases, i.e., the spectrum 
becomes softer An intiiguing case is V0332+53 (pink squares), 
whose behaviour differs from the other sources by displaying a 
positive correlation in three out of the four panels of Fig. This 
peculiar behaviour of the colours of V0332+53 was already re- 
ported in Paper II. Clearly, the energy spectra of BeXBs above 
10 keV are more complex than a simple power law. 

4.2. Spectral variability 

Although the energy spectra of all sources can be fitted with the 
models outlined in Sect. 13.21 the different components (power- 
law, cutoff, cyclotron line) fit different parts of the spectrum 
in different sources. Hence a direct compaiison of the actual 
value of some of the spectral parameters (i.e., photon index) 
between different sources may be misleading. Nevertheless, be- 
cause the spectral parameters vary smoothly along the outbursts. 
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Fig. 9. Variation of the energy of the cyclotron line with lumi- 
nosity (3-30 keV). Filled circles correspond to the HB and open 
circles to the DB. In 4U 01 15+63, red squares represent the re- 
sults from the "bump" model, whereas circles are the best-fit cy- 
clotron energy with the POWERLAWxHlGHECUT description 
of the spectral continuum. 

the search for correlations among different parameters and also 
with X-ray flux in each individual source is meaningful. Figure|5] 
shows some representative spectra for various flux levels of the 
HB and DB. The quoted X-ray luminosity in this figure was de- 
rived for the 2-100 keV range and normalised to the Eddington 
luminosity for a 1 .4 Mq neutron star. See Table[T]for the assumed 
distance estimation. 

The following general results are common to all or most of 
the sources: 

- The photon index anticorrelates with X-ray flux in the HB 
and correlates with it in the DB (Fig. |6]l. This means that as 
the flux increases, the spectrum becomes harder in the HB 
and softer in the DB. In either branch, these coiTelations are 
strong and significant. The Pearson's coiTelation coefficient 
for the variation of the photon index with flux range is p>0.9 
in most cases. Although correlations between the power- 
law photon index and the X-ray luminosi ty have been re- 
ported in the past (iReynolds. et al.l 1993; De vasiaetal .120111; 
iKlochkov et alJl2011h . this is the first time that a change in 
the slope of the correlation is seen, in coiTespondence with 
the position of the source in the HID. For 4U 01 15-1-63, the 
results from the two different ways to describe the spec- 
tral continuum are plotted in Fig. |6] Circles correspond to 
the POWERLAWxHlGHECUT model and squares to the 
CUTOFFPL+GAUSS (or bump) model, where GAUSS is a 
broad emissi on Gaussian profile that accounts for the "10- 
keV feature" ( Muller et al.ll20T2l) . 

- The cutoff energy approximately follows the same trend as 
the photon index, namely, it increases as the flux decreases 
in the HB and increases as the flux increases in the DB. As a 
result there is a positive correlation between the photon index 
and the high-energy cutoff: softer spectra correspond to high 
values of the cutoff energy (Fig.|7]i. 

- The fluorescent iron Ka line feature that results from repro- 
cessing of the hard X-ray continuum in relatively cool mat- 
ter is ubiquitous in all sources. Near neutral iron generates 
a line centered at 6.4 keV. As the ionisation stage increases 
so does the energy of the line. However, the energy separa- 
tion is so small that even ionised iron up to Fe XVIII can 
be thought as part of the 6.4 keV blend dEbisawa et al ]fT996t 
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Fig. 10. Representative power spectra of the horizontal and di- 
agonal branches at various flux levels for the fast-rotating pul- 
sar 4U 0115+63 (up), EXO 2030+375 (middle) and the slow 
pulsator 1 A 0535+262 (down). Note the shift toward higher fre- 
quencies at high flux. In lA 0535+262, Lq is not seen at low 
luminosities because it lies outside of the frequency interval con- 
sidered. Other features such as Llf in 4U 01 15+63 or a QPO at 
~40 mHz in lA 0535+262 are shown. The Lorentzian profiles 
of the pulse peaks and the QPO were removed for clarity. 



lLiedahlll2005h . The central energy of this component does 
not vary significantly during the outbursts and it is consis- 
tent with "near-neutral" material, that is, with the composite 
of lines from Fe II-Fe XVIII near 6.4 keV. The flux of the 
line increases with the continuum X-ray flux, indicating that 
as the illumination of the cool matter responsible for the line 
emission increases, so does the strength of the line (Fig.|8). 
In general, the line equivalent width remained insensitive to 
luminosity changes. The strong correlations of Fig. |8]imply 
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that the contributio n from thermal hot pla sma located along 
the Galactic plane (lYamauchi et al.ll200 9') is not significant. 
We would expect the ridge emission to affect the data only 
for very f aint observations, i. e., at the start and end of the 
outbursts dEbisavya et al 120081) . The small deviation from the 
linear trend with flux at very low fluxes would be consistent 
with this statement (Fig.[8]l. 

There is no universal trend between the energy of the 
cyclotron line and X-ray luminosity (Fig. |9)- Previous 
studies h ave show n a strong anticorrelation in V0332+53 
^Tsvgank ov et al.| [2010), someho w weaker anticorrelatio n in 
4U O'l 15+63 (i Miharaet al.l l2004: Nakajima et al." 2006|) and 
no con-elation in lA 0535+262 (Caballero e t al 2Qq3 ). We 
confirm the decrease of the cyclotron line energy with lumi- 
nosity in V0332+53, but we do not find a smooth anticorrela- 
tion in 4U 0115+63. When the POWERLAWxHIGHECUT 
model is used, the energy of the cyclotron line remains fairly 
constant within each branch, but its value is higher in the HB 
(~ 16 keV) t han in the DB (~ 11 ke V). T his resul t agree s 
with that of 'Tsvg ankov etal] ( l2007h and iLi et all (1201 2h . 
who reported a discontinuity in the E^yc - Lx relationship 
of this source with larger values of the cyclotron line en- 
ergy at low flux. However, if the continuum is fitted with the 
CUTOFFPL model and a broad emi ssion Gaussian, then no 
correlation is seen, in agreement with Miiller et al 1 (l2012h .In 
Fig- HI red squares represent the best-fit values of the energy 
of the cyclotron line using the "bump" model, while circles 
represent those obtained with the high-energy exponential 
cutoff model. 

Our results also show no clear trend between the en- 
ergy of the cyclotron line and luminosity in lA 1118-616, 
EXO 2030+375, and XTE J0658-073. However, we report, 
for the first time, a positive correlation between the cy- 
clotron energy and X-ray luminosity in lA 0535+262 on 
long timescales (Fig. |9l). A correlation h ad been reporte d 
only in the pulse-to-pulse analysis (.Klochkov et alj [201 lb . 
Previous studies did not find any eviden ce for variability 
dTerada et al.1 l2006t ICabaUero etaTI l2007h . However, note 
that these works show observations of lA 0535+262 when 
the X-ray luminosity was <1 x 10^^ erg s"' {LxlL^id < 0.1). 
As can be seen in Fig. |9] the energy of the cyclotron line 
begins to increase above this valu^l We verified that this re- 
lationship cannot be due to any artificial coiTelation between 
the line energy and the other parameters of the model by plot- 
ting the contours of the confidence regions of E^^y^ and any 
other parameter. Additionally, a linear correlation analysis 
gives a linear correlation coefficient of 0.89 and a probabil- 
ity that the points are not correlated of 6 x 10"**. 
To ensure that the Lorentzian profile was fitting the CRSFs 
and not the continuum, we carefully verified that the cor- 
responding width had reasonably narrow values. In all the 
cases, we obtained width values compatible with previ- 
ous works. For EXO 2030+375 and 4U 0115+63, which 
displays the lowest ener gy CRSF, the obta i ned width is 
of di e order of 3 keV (IWilson et al] 120081: iMiiUer et all 
l20T2h : for V0332+53 and lA 05 35+262, the retrieved 
width was of the ord er of 6-8 keV (iTsveankov et al.l[2010t 
ICaballero et al.' 2007); lA 1118-616 and XTE J0658-073, 
the systems that display the highest energy CRSF, also show 



^ In fact, the two data points abo ve ~ 1 x 10^^ erg s ' in Fig. 4 of 
iTerada et all ( l2006h and in Fig. 5 of ICaballero et al.l (120071) do not rule 
out an increase of the cyclotron line energy at higher luminosity. 



the largest width, o f 12-15 keV (iDoroshenko et al.1 120I0I: 
iMcBride at al.ll2006h . 
- The X-ray spectrum of the accreting pulsars analysed in this 
work is characterised by a power law (F = - 1) with an 
exponential cutoff in the range 10-20 keV. However, the 
spectral continuum may be very distorted by the presence 
of cyclotron resonant scattering features. Unlike BHBs or 
anomalous X-ray pulsars where the presence of different 
components (disc blackbody and power law, or two power 
laws) causes a sudden break in the continuum, the X-ray 
spectral continuum of accreting pulsars does not show any 
evidence for an abrupt change below or above 10-20 keV 
(Fi g. HI), despite the fact th at p hysical models such as those 
by iBecker & WolfH (l2007h and iFerrigno et al.1 ([2009) show 
that several well-defined physical components {i.e., thermal 
and bulk Comptonisation, blackbody emission at the base of 
the accretion column) are present. 

4.3. Timing aperiodic variability 

We conducted a systematic analysis of the rapid aperiodic vari- 
ability of the sources. Our aim is to provide a consistent descrip- 
tion of the power spectra of all sources, investigate whether the 
changes in spectral state are also evident in the power spectra, 
and study the evolution of the noise parameters throughout the 
outburst. 

We consistently used the same model to fit the power spectra 
of all sources. We found that the majority of power spectra are 
well represented by the sum of up to three Lorentzian profiles 
that were termed as L,, where / = 0, 1,2. Some power spectra 
displayed more complex substructure, which mainly consisted of 
the presence of wiggles without any preferred frequency range 
or narrow, generally weak, features confined to only one bin. 
Although these features may lead to a formally unacceptable 
fit (x^ > 2) in some cases, no attempt to coiTect for this ef- 
fect was canied out because it does not affect our results sig- 
nificantly. Here we focus on the broad-band noise components, 
that is, power covering relatively large frequency intervals. The 
use of Lorentzian profiles is particularly suitable for this kind 
of studies as it has been demonstrated in b lack-hole binaries 
(iPottschmidt et alJl2003l ; lAxelsson et al.ll2005h . Although it can- 
not model very narrow features, it provides a simple but con- 
sistent way to fit the main characteristics of the power spectra 
across different states and to track changes occurring on short 
timescales. In addition to the broad-band noise represented by 
Li, the fast rotating pulsars 4U 0115+63 and V0332+53 show 
other narrower components, whose characteristic frequency does 
not vary with luminosity. For a more detailed analysis of these 
components the reader is referred to Paper II. The peaks that 
correspond to the spin period and its harmonics were also fitted 
with Lorentzian profiles but the centroid frequency and width 
were fixed. Figure [10] shows some representative power spectra 
at various flux levels. 

Lo is a zero-centred Lorentzian that accounts for the noise 
below 0.05 Hz. Li is the main noise component and accounts 
for the noise in the range 0.1-1 Hz. It is also a zero-centred 
Lorentzian whose characteristic frequency increases as the flux 
increases. Its nns is always larger than 15%. L2 accounts for 
the broad-band noise at higher frequencies. Its characteristic fre- 
quency peaks typically in the range 1-5 Hz, while the nns is nor- 
mally below 15%. L2 appears at high flux, i.e., near the peak of 
the outbursts. This component is normally not statistically sig- 
nificant at the end of the outbursts. Although in most cases it 
shows up as a zero-centred Lorenztian, it may turn into peaked 
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noise, especially at the highest flux. Because of the narrower fre- 
quency range covered by the Lq component and typically larger 
error bars of low-frequency points, the addition of Lq, especially 
for the first appearance of this component, did not generally im- 
proved significantly. Hence, the addition of Lo was based on 
visual inspection of the residuals. In contrast, the introduction 
or removal of L2 was based on whether its presence or absence 
resulted in a significant improvement of the best fit;^^^. 

The following general results are common to all or most of 
the sources: 

- The characteristic frequency increases as the X-ray luminos- 
ity increases (Fig.fTTTi. 

- Whereas the overall 0.01-10 Hz rms does not change sig- 
nificantly within a branch, the fractional amplitude of vari- 
ability of the main component (Li) tends to decrease as the 
X-ray luminosity increases (Fig. [12}. With the exception of 
EXO 2030-1-375, this trend is more distinct for sources ex- 
hibiting two branches. That is, the source tends to be more 
variable in the HB. 

- L() and L2 are usually not present in the HB. The disap- 
pearance of Lq in the HB is probably due to the fact that 
this component falls below the frequency interval considered 
(vo < 0.01 Hz) at very low fluxes in agreement with the ob- 
served frequency shift (Fig.fTTTi. 

- Whenever Li and L2 appear simultaneously in the power 
spectrum, their frequencies correlate (Fig. [13). 

The saturation at high luminosity seen in the luminosity- 
frequency diagram of some sources (most notably in 



XTE J0658-073, Swift J1656.6-5156, and EXO 2030-H375) 
is due to the appearance of the L2 component. The correlation 
between the characteristic frequency and the flux is broken when 
this extra component is added because it prevents Li from shift- 
ing toward higher frequencies. As an exercise, we removed the 
L2 component from the power spectra of Swift J1656. 6-5156. 
The absence of L2 gives worse fits, but vi shifts up (squares in 
Fig.fTTI) and the correlation holds even at the highest luminosity. 

4.4. Correlation between spectral and timing parameters 

The spectral parameters correlate tightly with the X-ray flux. In 
some sources, though, this correlation changes sign at very low 
fluxes (in the HB). The characteristic frequency of the broad- 
band noise also shows a smooth relationship with flux, albeit 
with more scattering. The question is then whether the spectral 
and timing parameters correlate with each other. Figure [l4ldis- 
plays the Li maximum frequency as a function of the photon in- 
dex. These two parameters seem to correlate in most sources, al- 
though the sign and strength of the correlation vary significantly. 
We shall discuss the implication of this result in Sect. 15.31 

5. Discussion 

The main goal of this study is to perform a detailed spectral and 
timing analysis of accretion-powered pulsars with Be compan- 
ions in an attempt to characterise this type of systems as a group. 
More precisely, we wish to investigate i) whether accreting X- 
ray pulsars display spectral states as black-hole and low-mass 
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X-ray binaries do and characterise those states, and ii) to search 
for coiTelations between spectral and timing parameters during 
major outbursts in an attempt to constrain the accretion models 
in X-ray pulsars. 

In this section we discuss the implications of our work. We 
interpret our results in the context of the model that suggests 
the existence of tw o different accreti on regimes, defined by a 
critical luminosity (iBecker et al.ll201Z) . First, we make a short 
introduction of the spectral formation mechanism in accreting 
pulsars and then we summarise and discuss our results. 

5.1. Spectral formation 

A detailed description of the emission properties of an accretion- 
powered pulsar is a complex task. First, it requires an under- 
standing of the processes by which mass is transferred, captured, 
deposited on the neutron star surface in the form of an accretion 
column and converted into high-energy radiation. Second, this 
radiation has to find its way out, which involves a knowledge of 
the interactions of the X-rays produced close to the neutron star 
surface with the highly magnetized plasma forming the magne- 
tosphere. Third, the X-rays emanating from the polar caps are 
subjected to absorption and reflection processes due to the am- 
bient matter. 

This complexity translates into a lack of a fundamental phys- 
ical model that yields results that agree with the observations. 
The common practice is then to fit phenomenological multi- 
component models to the energy spectra of accretion-powered 
pulsars. The parameters of these model components are diffi- 
cult to interpret physically. Some attempts to allevi a te th is sit- 
uation have been made by iBecker & V^ll (I2005L I2007L see 
also Ferrigno et al. 2009), who developed a new model for the 
spectral formation process in X-ray pulsars based on the bulk 
and thermal Comptonisation of photons due to collisions with 
the shocked gas in the accretion column. The accretion flow is 
channelled by the strong magnetic field into the polar caps, cre- 
ating an accretion column. Most of the photons are produced 
at the base of the column, just above the neutron star surface 
(thermal mound). The low-energy (blackbody) photons created 
in the mound are upscattered due to collisions with electrons 
that are infalling at high speed (bulk Comptonisation). Thermal 
Comptonisation (in this case high-energy photons lose energy) 
also plays a significant role in shaping the spectrum and it is 
responsible for the formation of the exponential cutoff. 

5.1.1. Two accretion regimes 

There exists a critical luminosity, L^nt , at which the deceler- 
ation of the accreting flow to rest at the neutron star surface 
changes from being do minated by radiation pressu r e to occur via 
Coulomb interactions (iBasko & SunvaevI Il976h . IBecker et al.l 
suggested that these two different regimes can explain the 
bimodal behaviour of the variability of the cyclotron energy with 
luminosity. Some sources, such as V0332-I-53 display a negative 
correlation between the cyclotron energy and the source lumi- 
nosity, others, such as Her X-1 (Staubert et al. 2007), show a 
positive correlation. 

Two accre t ion r egimes have also been invoked by 
iKlochkov et alj (1201 lb to explain the pulse-to-pulse variability 
in some X-ray pulsars. It is illustrative to compare our results 
with those of Klochkov et al. (2011) because of the different 
timescales involved. We have studied the long-term variability 
patterns of accreting pulsars by sampling timescales of X-ray 




Fig. 12. Variation of the fractional amplitude of variability of 
component Li with luminosity. Filled black circles correspond 
to the HB and open circles to the DB. 
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variability of the order of 1000 s. iKlochkov et al.l (1201 lb stud- 
ied the pulse-to-pulse spectral variability, i.e., timescales of the 
order of the pulse period (f spin ~ 1 - 100 s). Their sample 
included four X-ray pulsars, three of which appear in our list: 
V0332-H53, 4U 01 15+63, and 1 A 0535+262. To achieve enough 
photon statistics they restricted their analysis to high-flux states. 
In the context of the present work, this means that they analysed 
data of the DB in V0332+53 and 4U 01 15+63 and of the high- 
flux part of the HB in 1 A 0535+262. Despite the order of mag- 
nitude difference in the variability range, the same correlations 
are observed. They found that the photon index of the power- 
law component correlates with X-ray intensity in V0332+53 and 
4U 0115+63, but anticorrelates in lA 0535+262, which is the 
same result as our Fig|6] Likewise, they reported a coiTelation 
between the energy of the fundamental cyclotron line with count 
rate in 1 A 0535+262, as we do with our long-term pulse-average 
spectral analysis (Fig.|9]l- 

5.2. Source states 

This work shows that, as BHBs and LMXBs, accreting X-ray 
binaries with high-mass companions also exhibit spectral states 
in the hardness-intensity diagrams. At low and intermediate 
X-ray luminosities the sources populate the horizontal branch 
(HB), whereas at higher luminosities the sources trace a diag- 
onal branch (DB). However, unlike BHBs and LMXBs, where 
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the hardness of the source is a distinguishing property of the 
accretion states (the so-called hard and soft states), the spectral 
branches in accretion-powered pulsars roughly sample the same 
colour interval. If a source displays the two branches, then the 
timescales of the colour and flux changes are faster in the HB 
than in the DB. The source covers the HB in hours to few days 
and the DB in days to weeks. In contrast, systems that show only 
the HB can remain in that branch for weeks. 



With flie exception of 4U 0115H-63 and V0332-H53 (see 
Paper II), there is no clear evidence for hysteresis in the HID, 
i.e., the same colour corresponds to two different values of the 
count rate depending on whether the source is in the rise or 
decay. This lack of hysteresis could be due to the fact that, in 
most cases, the rise of the outburst is much badly sampled than 
flie decay (see Fig.[B. For Swift J1656.6-5156, lA 0535-1-262, 
lA 1118-616, and GRO J1008-57, no or very few data points 
of the rise were observed. On the other hand, in KS 1947-1-300, 
which is the source with the most complete coverage of the rise, 
no indication of hysteresis is found. To check this result, we per- 
formed a colour analysis on the ASM light curves and we did 
not find significant deviations between the rising and decaying 
tracks. However, it should be noticed that the errors of the ASM 
colours are, in the best case, roughly of the same size as the am- 
plitude of the hysteresis effect. It is also worth mentioning that 
while we did not observe significant hysteresis in the HID of 
EXO 2030-H375, the characteristic frequency of the broad-band 
noise in this source exhibits lower values during the rise than 
during the decay (Fig.fTTTi. 



One of the main findings of this work is that the transition 
from the HB to the DB occurs when the source luminosity in- 
creases above a certain value, which we find to be in the range 
~ 0.06 - 0.3LEdd (~ 1 - 4 X 10" erg s"'). We observe that this 
value is different for different sources. We suggest that the lu- 
minosity at which the transition occurs is related to the criti- 
cal luminosity and that each spectral branch can be associated 
with one of the two modes of accretion proposed by Becker et al] 
120121) . The DB, which corresponds to a high-luminosity state, 
would be associated with the supercritical mode and the HB with 
the subcritical one. 

The case of KS 1947-1-300 and lA 0535-H262 would cor- 
roborate this association. These two sources reached similar 
peak outburst luminosity. However, only KS 1947-t-300 shows 
the two branches. KS 1947 -h300 reached a peak luminosity of 
0.4LEdd and made the transition from the HB to the DB at 
~ 0.06LEdd o r ~ 1 X 10^^ er g s"', assuming a distance of ~ 10 
kpc (Neguer uela et al.ll2003i) . In contrast, lA 0535H-262, with 
the same peak luminosity stayed in the HB during the entire 
outburst. We conclude that the intensity of the source is not the 
only parameter triggering the change of state and some other pa- 
rameter must play a role too. The magnetic field appears as the 
best ca ndidate to drive the transition. According to Becker et al] 
(1201 2h . the critical luminosity Lent depends on the magnetic 
field, or equivalently, on the energy of the cyclotron line. The 
higher the magnetic field, the higher t he critical lumino sity. For 
typical neutron star parameters (see Beck er et al.ll2012l for de- 
tails), the critical luminosity can be estimated as 



1.28 X 10"(£cyc/10keV)"'/'" erg s 



16/15 



(3) 
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where Ecyc is the energy of the cyclotron line. lA 0535+262 
displays one of the highest energy cyclotron lines, at ~ 50 keV 
(Caballero et al. 2008). Substituting in the above equation yields 
La-it ~ 7 X 10^^^ erg s"', which is slightly larger than the out- 
burst peak luminosity. No cyclotron line has been reported for 
KS 1947+300, which would indicate a lower magnetic field, in 
agreement with the lower luminosity at which the spectral tran- 
sition occurs. The luminosity at which the transition is observed 
in KS 1947+300 would imply a cyclotron line at around 8 keV. 
Unless the cyclotron line is very strong such a feature would be 
difficult to detect because it would lie in a region where other 
spectral features dominate. First, the fluorescence iron line at 
6.4 keV and possible its associated absorption edge. Second, 
some X-ray pulsars show si gnificant wiggle re siduals at 8-12 
keV, whose origin is unclear (ICoburn et al.ll2002 ). Nevertheless, 
it is worth noticing that in the high-flux average spectrum of 
KS 1947+300, an additional absorption component at about 9 
keV provides an improved fit with respect to the model used for 
the individual spectra. The;\f^ decreases from 104 for 72 degrees 
of freedom to 74 for 69 degrees of freedom. This component can 
be modelled with an absorption Gaussian profile (GABS) and 
could be associated to a cyclotron line that in the single spectra 
would be too weak to be unveiled. The absorption could also be 
accounted for with an edge {x^-l\ for 70 degrees of freedom), 
although in this case the inspection of residuals reveals that this 
component is unable to properly fit the entire spectral region of 
interest. In short, it is difficult to tell whether the source actually 
displays a cyclotron line at this energy, but if so, its low charac- 
teristic energy would agree with the critical luminosity reached 
by the system. 

The larger magnetic field in lA 0535+262 might also ex- 
plain the reason that previous studies of this source did not find 
a correlation between the energy of the cyclotron line and the 
X-ray luminosity. Teradaetal. (2006) speculated that this en- 
ergy might begin to change at a higher luminosity if the object 
has a higher surface magnetic field. Here we report observa- 
tion s of 1 A 0535+262 at higher luminosities than previous stud- 
ies (Kendziorra et al."1 994t IWilson et alj|200l ICaballero et alj 
[2007; Terada et al. 200q). As it can be seen in Fig.|9] the energy 
of the CRSF only begins to increase significantly when the X-ray 
luminosity is above ~ 3 x lO-'^ erg s"'. 

XTE J0658-073 is another source whose peak luminosity is 
close to its critical luminosity (see Table[TJ. Although it does not 
display a proper DB, an attempt to move to that branch can be 
observed in Figs. |3]and|6] In particular, note the upward turn of 
the photon index at Lx ~ O.lLEdd in Fig.|6l 

iBecker et al.1 (l2012h also showed that as the luminosity in- 
creases, the height of the emission zone inside the accretion col- 
umn approaches the neutron star surface in the subcritical accre- 
tion state, while this emission zones moves upward in the su- 
percritical accretion state. Comptonisation occurs in the region 
between the radiative shock and the neutron star surface (the 
sinking region). In the supercritical state this region is relatively 
small , just a few tens of meters (see eq. (40) in iBecker et al.l 
120 12h . but its height increases as the luminosity increases. In 
this region, the effective velocity of the comptonising electrons 
is highly reduced because diffusion (outwards) and advection 
(inwards) are almost balanced. Hence photons would not ac- 
quire enough energy through bulk Comptonisation to populate 
the higher energy band and a softening of the spectrum would 
be expected with increasing luminosity. This agrees with the ob- 
servations: in the DB, the spectrum becomes softer (soft colour 
decreases, F increases) as the luminosity increases (Figs.|3]and 
|6]). In the subcritical mode, the typical emission height is a few 



kilometers (see eq. (51) in IBecker et al.ll20T2h but decreases as 
the X-ray luminosity increases. As the size of the sinking re- 
gion decreases with increasing luminosity, the optical depth in- 
creases, which results in harder photons. In the HB, the spectrum 
becomes harder as the luminosity increases. 

Finally, it is worth noticing that the four sources that dis- 
play the two states have short orbital (Poib < 50 days) and spin 
(^'spin << 100 s) periods, whereas the sources that do not tran- 
sit to the DB have, in general, PqA > 100 days or f spin ^100 
s (see Table [TJ. The result that systems harbouring fast-rotating 
neutron stars in narrow orbit s display tw o branches in the HID 
can be linked to the work by Reig (2007), who found that BeXB 
with short spin and orbital periods are more variable in the X-ray 
band, i.e., they exhibit higher amplitude of variability as mea- 
sured by the root-mean-square rms over the long-term (years) 
X-ray light curves. Reig (2007) explaine d this result in the con- 
text of the viscous decretion disc model (lOkazak i & Neguemelal 
2001). This model predicts the truncation of the Be star's disc 
due to the tidal interaction exerted by the neutron star. In the 
truncated disc model, the material lost from the Be star accu- 
mulates and the dis c becomes denser m ore rapidly than around 
an isolated Be star (lOkazaki et al.ll20()2h . Truncation is favoured 
in systems with short orbital periods and low eccentricities. In 
such systems mass transfer would occur only when the disc is 
strongly disturbed. Accretion of large amount of material from 
the distorted disc would give rise to very bright X-ray outbursts. 

5.3. Correlation of spectral and timing parameters 

We searched for correlations between the different model pa- 
rameters that describe the energy and power spectra. Two ap- 
proaches were taken to investigate the relation between the X- 
ray emission and the rapid aperiodic variability. In the first one, 
we tracked the spectral and timing parameters following an X- 
ray outburst, in which case the general trend of the relation with 
luminosity was inferred (Figs. |6] |8] |9] [TT] and [12). In the sec- 
ond approach, we directly compared parameters between them- 
selves (Figs.l7]andfl4li. To quantify and assess the statistical sig- 
nificance of the various relationships we obtained the Pearson's 
correlation coefficient. The results are summarised in Table |4] 

Of particular interest are the correlations between spec- 
tral and timing para meters. According to the model by 
IBecker & WolfiTJIOOTl) . the hard power-law component is due 
to bulk inverse Comptonisation in the accretion column. Low- 
energy photons are upscattered in the shock and eventually dif- 
fuse through the walls of the column. Therefore, the photons that 
populate the hard power-law tail come from a region close to the 
surface of the neutron star (length scales ~ 10^ cm). 

On the other hand, aperiodic variability must be generated at 
much larger distance. The evidence in support of this statement 
comes from the correlation between the characteristic frequency 
of the broad-band noise and the X-ray flux and from the milli- 
Hertz range of frequencies of the QPO seen in accreting X-ray 
pulsai's. 

The radius of the magnetosphere, r„„ depends on the mass 
accretion rate (iDavies&Pringlell 19811) 

r,„ = 1.2xl0XfV3i'(^) cm (4) 

where M15 is the mass accretion rate in units of 10'^ g s"' and 
the dipole magnetic moment in units of 10""' G cm"*. As the X- 
ray flux increases, presumably as a consequence of an increase 
in the mass accretion rate, the magnetospheric radius decreases. 
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If the processes that originate the broad-band noise are Unked 
to motion of matter outside the magnetosphere, then we would 
expect that as the magnetosphere shrinks, the characteristic fre- 
quency increases, as observed (Fig.fTTTi. Thus this correlation is 
consistent with and extra-magnetospheric origin of the aperiodic 
variability. 

Likewise, QPOs in HMXB lie typically in the milliHertz 
range (jJames et al. 2010). If one assumes that QPOs are pro- 
duced as a result of Keplerian motion of inhomogeneities in 
an accretion disc, this frequency range (0.01-1 Hz) agrees with 
the QPO being originated outside (but near) the magnetosphere, 
i.e., at length scales r,„ > 10'^ cm. Indeed, a Keplerian fre- 
quency of vqpo = 100 mHz would correspond to a radius 
r = (GM/47T^vlpgy'^ ~ 7.8 x 10** cm. Note also that in LMXBs, 
peaked (broad-band) noise have been seen to developed into 
QPOs, indicating a physical relationship between broad-band 
noise and QPOs (van der Klis 2006). Thus it is reasonable to as- 
sume then that the aperiodic variability originates from the same 
physical region as the QPOs. 

Since the power-law photons and those responsible for the 
aperiodic variability originate in very different physical regions, 
we would not expect the spectral and timing parameters to cor- 
relate. We have plotted the relationship between the power-law 
photon index and the frequency of the main noise component 
Li in Fig. [141 Despite the different location of the production 
mechanisms, a distinct correlation between the power-law index 
and characteristic frequency is seen in almost half of the sources. 
Mimicking the transition from a negative to a positive correlation 
seen in the Lx -T diagram as the flux increases (Fig. |6]l, the sign 
of the correlation between vi and F also changes, from negative 
in the HB to positive in the DB (Fig. [T4b . Whether this result 
is simply the consequence of the strong correlation between the 
characteristic frequency and the X-ray luminosity (Fig.lTTI). i.e., 
the frequency acts as a proxy for luminosity, or else it has a more 
profound physical meaning remains to be solved by further stud- 
ies. 

6. Conclusions 

We have investigated the correlated X-ray timing and spectral 
variability of accreting X-ray pulsars with Be-type compan- 
ions during major X-ray outbursts. The aim was to investigate 
whether this type of systems exhibit source states and in case 
they do whether the X-ray timing and spectral properties are 
strongly correlated as it is seen in neutron-star low-mass and 
black-hole X-ray binaries. Although there may be some dis- 
crepant behaviour in particular sources, which substantiates the 
complexity of the accretion process, we were able to extract gen- 
eral patterns of variability. 

The evolution through the hardness-intensity diagram sug- 
gests that in the early and late phases of the outburst, BeXBs un- 
dergo state transitions. As the source evolves along the outbursts 
it transits from the horizontal branch to the diagonal branch and 
back. We show that the state transitions occur when a critical 
luminosity is reached. Only sources whose peak luminosity is 
well above the critical limit do exhibit the two branches. Lower 
than the critical luminosity sources display only the horizontal 
branch. Because the value of the critical luminosity depends on 
the pulsar magnetic field, the luminosity at which the transition 
takes place varies across different sources. For typical values of 
the magnetic field in Be/X-ray pulsars, the critical luminosity is 
ex pected to be of the ord er of a few times 10^^ erg s"' . 

iBecker et al. I (I20TI showed that two different accretion 
regimes can explain the bimodal behaviour of the cyclotron line 



with luminosity displayed by diff'erent sources. We propose that 
the two branches correspond to these two accretion regimes, 
likely to be related to the way in which the accretion flow is 
decelerated in the accretion column: radiation pressure in the su- 
percritical regime (diagonal branch) and coulomb interactions in 
the subcritical regime (horizontal branch). In this work we show 
that the continuum traces the two regimes as well. Despite the 
complex spectral continuum, the power-law photon index cor- 
relates with X-ray luminosity and also exhibits two branches. 
When plotted as a function of X-ray luminosity, the correlation 
changes sign in coincidence with the change of state. In the hori- 
zontal branch, the power-law photon index decreases with X-ray 
flux, while it increases in the diagonal branch. We speculate that 
this behaviour is due to the different dependence of the accretion 
column height with flux in the two accretion regimes. 

In contrast, although a clear positive correlation between the 
frequency of the broad-band noise components with flux is ob- 
served, the frequency does not trace the two branches. This result 
would support the idea that the aperiodic variability originates 
further away from the neutron star surface, outside the accretion 
column. 

We report for the first time a positive correlation be- 
tween the energy of the cyclotron resonant scattering feature in 
lA 0535+262 with luminosity on long time scales. The energy 
of the cyclotron line is seen to start to increase at a higher lumi- 
nosity than other sources, presumably due to its larger magnetic 
field. This is only the second X-ray pulsar (after Her X-1) to 
display such positive correlation. 

Finally, some sources show a correlation between the pho- 
ton index and the characteristic frequency of the aperiodic noise 
components, implying that the accretion column, where energy 
spectra are generated somehow communicates with the inner ac- 
cretion disc, where the aperiodic variability is supposed to orig- 
inate. This result imposes a tight constrain to the models that 
seek to explain the spectral and timing variability in accretion- 
powered pulsars. 
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